Electrically tunable two-channel Kondo fixed points in helical liquids 
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We study a quantum dot coupled to two edge states of a quantum spin Hall insulator through 
electron tunnelings in the presence of a Rashba spin-orbital interaction induced by an external 
electric field. We show that if the electron interactions on the edge states are repulsive, there are 
two possible phases, depending on the Luttinger liquid parameter K. For 1/2 < K < 1, the low- 
temperature physics is controlled by a previously identified two-channel Kondo fixed point. For 
the edge states with even stronger repulsive interactions, i.e. 1/4 < K < 1/2, the system reaches 
another phase at low temperatures, described by a new two-channel Kondo fixed point. This phase 
is separated from the original one by a continuous phase transition upon varying the value of K 
through the external electric field. The corresponding quantum critical point is described by a free 
Dirac fermion backscattered by a local potential. We investigate the low-temperature properties 
associated with this new fixed point and also discuss the scaling behaviors of the system in the 
quantum critical regime. 



PACS numbers: 71.10.Pm 72.10.Fk 72.70.-|-m 

Introduction. — The experimental discoverjsi of the 
quantum spin Hall (QSH) insulator in HgTe quantum 
wells following its theoretical predictions^ paved a new 
road in the study of topological phases. The QSH stated 
is a member of the topologically nontrivial states of 
matter with the symmetry protected topological order 
These states of matter have a finite bulk gap and, in the 
mean time, support gapless edge or surface excitations. 
For the QSH insulator, the edge states are helical in the 
sense that, on each edge, there is a counter-propagating 
Kramers' pair of states with opposite spin polarizations. 
The stability of the helical edge states against poten- 
tial scattering is protected by the time-reversal (TR) 
symmetryi^^— Contrary to the case of simple potential 
scattering, the edge electrons can backscatter from the 
magnetic impurities through spin exchange, and thus the 
TR symmetry can no longer protect the helical states 
from mixing. The one-channel Kondo (ICK) effect of 
helical edge states was studied in Ref. I9l-[l2i 

In a recent work, a quantum dot (QD) coupled to 
two helical edge states was studiedji^ It is well known 
that such a system realizes the usual ICK effect for non- 
interacting electrons on the edgesJ^ii^ In Ref. |T^, it 
is shown that for weakly repulsive interacting electrons 
on the edges, with the Luttinger liquid (LL) parameter 
K < 1, the system is driven to a two-channel Kondo 
(2CK) fixed point. This result is non-trivial since in the 
context of a QD coupled to two LL leads, a much stronger 
Coulomb repulsion (K < 1/2) is warranted in order to 
realize the 2CK physics. 

A more complete description of the 2CK physics in 
a QD coupled to two helical edge states must take 
into account the Rashba spin-orbital interaction because 
this interaction, which can be tuned by an external 
gate voltage, is a built-in feature of a quantum wellJ^ 



Moreover, the HgTe quantum wells exhibit some of the 
largest known Rashba couplings of any semiconductor 
hetcrostructuresji^ In fact, it was found very reccntlj*i^ 
that the presence of a Rashba coupling has profound ef- 
fects on both the Kondo temperature and the transport 
properties of the helical liquids in the presence of a single 
magnetic impurity. 

In this work, we study the effect of the Rashba cou- 
pling on the system studied in Ref. [ij. We find that if 
the electron interactions on the edge states are repul- 
sive, there are two possible phases, depending on the 
value of K. For 1/2 < K < 1, the low-temperature 
physics is controlled by the 2CK fixed point identified 
in Ref. [Tsl . In this region, we find that the impurity 
entropy at zero temperature Simp and the temperature 
dependence of the tunneling conductance between the 
two edges arc not affected by the presence of the Rashba 
coupling. On the other hand, the Rashba coupling re- 
veals its presence through the temperature dependence 
of the impurity specific heat Cimp and the dynamical 
structure factor Si{uj) of the spin in the dot at low tem- 
peratures, i.e. C^p cx T^+'^/K-2 foj. < K < 1, 
Crnp oc T4^-2 fo^. i/2<K < 1/^3, and S^{0) cx T^-^ 
with i — y, z. More interestingly, we show that the 
system reaches another phase at low temperatures, de- 
scribed by a new 2CK fixed point when 1/A < K < 1/2. 
In this region. Simp = ln(2/^), Cimp « 
5i(0) cx r2/f-i^ rpj^p phases are separated by a 
quantum critical point (QCP), which is described by a 
free Dirac fermion backscattered by a local potential. In 
the following, we outline the derivation of these results. 

Model. — We consider two helical edge states of a 
QSH insulator brought close to each other at a tunneling 
junction. A QD is then placed at the middle of the junc- 
tion and coupled to the two edge states through electron 
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tunnclings. When the number of electrons in the dot is 
odd and the repulsive interaction between the electrons 
in the dot is much larger than the tunneling amplitude, 
the system can be described by the Kondo Hamiltonian 
H = Ho + Hi + HKr^'^ where 

™ — 1 



= ($1 ± $2)/y2, and e,/a = (Gi ± es)/^ 
the Hamiltonian for the edge states can be written as 

Ho + H, =Ea=.,a f/'^^[^(5-0")' + where 
the LL parameter K and the speed of the collective ex- 
citation V depend on both the Coulomb interaction and 
the Rashba coupling strength. On the other hand, Hk 
becomes 



and 



m=1.2 
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In the above, m — 1,2 label the two edges, ^E'm = 

[V'm+jV'm-]*, V'm+ = V'mit' "^rn- = ^m/?i, Cr =t, |= 

, and 5 is the spin operator for the spin- 1/2 impu- 
rity. Here, Hq and Hi denote the kinetic energy (includ- 
ing the Rashba coupling with the strength a) and the 
electron Coulomb interaction of the helical edge states, 
respectively. Hk describes the Kondo interaction be- 
tween electrons in the edges and a spin- 1/2 magnetic im- 
purity at a; = 0. The magnetic anisotropy, i.e. Ji± ^ Jiz 
with / = 1, 2, is induced by spin-orbital couplingi^ 

The Rashba term in Hq can be absorbed into the ki- 
netic term by a spinor rotation = e~*2'^^\If,„^i^ 
where 9 = tan^^ (a/vp). By rotating also the 
impurity spin S ~ e*^'^^Se~*^'^^ , Hq bec omes 
'^Fj2m=i.2l'^^ *}„o'3«c^x*m , where VF = \/vp + a^, and 



H^ 



K 



m— 1,2 i—x,y,z 

+ E 

m=l,2 



(3) 



m^^n— 1,2 



where J;^ = J;j_, Jiy = Jji^ cos^ 6* + sin^ 6*, = 
Jiz cos^ 9 + Ji± siii^ 9, and Jie = {Ji± — J/z) sin cos 
with ; = 1,2. 

To proceed, we bosonize the Kondo Hamiltonian H = 
Hq + Hi + Hk according to the formula: 'ipm± ~ 
1 gT»\/47r0m± 22 ^i^QYe an IS the short-distance cut- 

\/27rao 

off. Since the Rashba coupling respects the TR sym- 
metry, we expect that the helical liquid on each edge 
is still described by the LL. Therefore, by defining the 



bosonic fields $r 



Hk = 



JlxSx 



■ COS 



TTOo 

JlySy + JieSz . 



27r$^(0) cos %/2^$a(0) 

2^$^(0)j cos[V27r$a(0) 



J2xSx 



Trao 

J2ySy + J2eSz 
TTOo 

J2zSz + J2ESy 



27r<I>s(0) COS V27rea(0) 



sm 



sm 



1 



TTflo 

{JlzSz + JlESy)dxQsiO) 





COS 




V^*a(0) 


sin 
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In general, there are extra local backscattering terms 
caused by the QD in the Hamiltonian H. However, these 
terms are irrelevant as long a.s K > 1/44 In the follow- 
ing, we shall focus on the regime with 1/4 < X < 1, so 
that these local backscattering terms can be neglected 
for low-energy physics. 

Scaling analysis. — Near the Gaussian fixed point 
{Jii = = JiE with I = 1,2 and i = x,y,z), the 
scaling dimensions of the various terms in Eq. Q are 
A(Ji,) = K = A(Jiy) = A(Jib), A(J)i, = 1, and 
A(J2,) = K/2 + 1/{2K) = A(J2y) = A( J2z) = A( J2£). 
Thus, for K < 1, the J2 terms decrease while the Ji 
terms grow when the temperature is lowered. 

In order to study the physics at the strong-coupling 
regime when Ji is of order 1, we follow Ref. [Tj by 
employing the Emery-Kivelson unitary transformation 

U = exp 
and Qs/a 



iV2TiK<l>,{0)Sz ^ where = <i>s/a/VK 



KG. 



la- 



The transformed Hamiltonian 



H = UHW is of the form 



— cos ^/2^TK^a{fi) 




■ cos 



^|e.(o) 



V27rA'$^(0) cos V2^$a(0) +5H,{b) 



where i/^^^^ = §/dx[(a,e,/,)2 -f (9,$,/ J^], \ = ^ 

\J2'kKv, Jl± = {Jlx + Jly)/2, SJl± = Jlx - Jly = (J/± 



3 



Ji,) si-a^e {I = 1,2), and 

6H = 51 1 5*:^ cos y/S-KK^siO) 



Sy sin 



V87rA'$^(0) I 



X cos 



V2^$a(0) 



92 COS 



|e.(o) 



X^Sx COS 

+53^2 sin 



V87rA"$^(0) -S'ysin \/8ttK^s{0) } 



\/27rA'$40) 



cos 



+54S'z sin 
+ff5|<5'y cos 



\/27rA'l>a(0) 
V27rA'$,(0) 



sm 




5:, sin V27rA'l>^(0) | 



xa,e,(0) +g6sin[V2^$a(0) 
x|5j;Cos V27rA:$,(0) 



5^ sin 



^|e.(o) 



\/2^$,(0) } 



27ra( 



,93 



Trao 



-,.94 = 



Trao ' 



In the above, gi = 52 
05 = and 06 = - — ■ 

At the vicinity of the point X = Jix = Jiy ~ Jiz = 
JiE ~ (/ = 1,2), we may calculate the scaling dimen- 
sions of the various J terms, and the results are shown in 
the first column in Tabic HI When 1/4 < AT < 1, among 
the possible relevant Ji±, J2±, Jie, and 171 terms, the Jij_ 
term is the most relevant one. At low temperatures, we 
expect that the effective Hamiltonian of H is described 
by the following fixed-point Hamiltonian 



= + H, 



(0) 



JlA 

Trao 



■ Sx cos 



V2^$,(0) 



(6) 



We notice that 77, is exactly the 2CK fixed-point 
Hamiltonian which has been thoroughly studied in Ref. 
[T3I . Since [Sx,H^] = 0, we may set Sx to be its eigen- 
values ±1/2, yielding 77± = 77,(5:, = ±1/2). In 7J±, 
the symmetric sector is described by a spinless LL, while 
the antisymmetric sector is described a spinless LL with 
an impurity backscattering term at x = and the ef- 
fective LL parameter K/2. The cosine term is relevant 
ion K < 1^ which cuts the antisymmetric sector into 
two separated pieces at x = at zero temperature. At 
finite temperature, tunneling between the two half wires 
is allowed, leading to the perturbations AiOi, where 
Oi = 4'^(0)\I'b(0) + H.C. and '^a{B) referred to fermions 
in the two separated pieces of the wire. The stability 
of this fixed point in the present case can be examined 
by calculating the scaling dimensions of the various per- 
turbations around it in a way similar to that given in 
Refs. [l^ and [23, and the results are shown in the second 
column in Table HI We see that for 1/2 < 7^ < 1 all per- 
turbations are irrelevant around 7J, , and this guarantees 
the stability of the 2CK fixed point in this region even in 
the presence of the Rashba interaction. 



Op.\FPs 


A — — Jlx(y) 


2CK 


2CK' 


Ji± 


K/2 


N/A 


N/A 


J2± 


1/(2K) 


1/A 


1/A 


JlE 


K 


A/2 + 1/(2A) 


5/(8A) 


A 


1 


1 + 1/(2A) 


1 + 5/(8A) 


91 


5A/2 


2A 


5/(8A) 


92 


27^' + 1/(2A') 


2A + 1/A 


1/A 


93 


A/2 + 1/(2A) 


A/2 + 3/(2A) 


13/(8A) 


94, 


A/2 + 1/(2A) 


3/(2A) 


13/(8A) 


95 


A/2 + 1 


A/2 + 1 


1 


96 


K + 1/(2/^) 


A/2 + 1/A 


1/A 


Ai 


N/A 


2/A 


2/A 


A2 


N/A 


N/A 


1/(2A) 



TABLE L Scaling dimensions of the operators (Op.) at dif- 
ferent fixed points (FPs), except those at the Gaussian fixed 
point, which are indicated in the text. 2CK refers to the 
fixed point with A = 0, |Jix| — !■ +c>o, and 2CK' refers to 
the fixed point with A = 0, | Jix |, liJJix | — )■ +00. J2,; = 
with i — x,y,z for all fixed points. Only the coefficients of 
corresponding operators are shown. 



On the other hand, for 1/4 < 7^ < 1/2, the 171 term be- 
comes relevant, and it renders 2CK fixed point unstable. 
Since the Jij_ term also fiows to strong coupling in this 
region, we expect that at low temperatures the system is 
described by the new fixed-point Hamiltonian 



= 7Jl")+^^5.cos 



+H, 



where = (cos 



27rao 

(0)+ -II 
Trao 



V87rA:$,(0) 



Sx cos 



\/27r7v$a(0) 



(7) 



We refer to this new 



V2^$a(0) 

fixed point as the 2CK' one in the following. Again, 
[Sx,H^] = 0, and we may set Sx = ±1/2 so that 
H± = Ht{Sx = ±1/2). Both the symmetric and the 
antisymmetric sectors in H± amount to a spinless LL 
with an impurity backscattering term at x = and the 
corresponding effective LL parameters Ks ~ 2K and 
Ka = K/2, respectively. For 1/4 < A' < 1/2, the J± and 
the 6Ji± terms both ffow to the strong-coupling regime 
and eventually cut both sectors into two separated pieces 
at a; = 0. At finite temperature, tunneling between the 
two half wires is allowed, leading to the perturbations 
AiOi and A2O2, where 62 = ^^(0)*b(0) + H.c. and 
^A(B) referred to fermions in the two separated pieces 
of the wire for the symmetric sector. The scaling dimen- 
sions of the various perturbations around the 2CK' fixed 
point is given in the third column in Table HI We see that 
all the perturbations are irrelevant for 1/4 < K < 1/2, 
except the term, which is marginal. This confirms that 
77, is indeed the low-energy effective Hamiltonian in this 
region. 
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Physical properties. — One of the most remarkable 
features of the multichannel Kendo effect is the existence 
of the fractionally degenerate ground state, which reveals 
itself in the impurity entropy at zero temperature. To 
compute it, we notice that Sx commutes with the fixed- 
point Hamiltonian, and thus we may write the partition 
function Z as Z = Z+ + Z- for 1/2 < K < 1 and Z = 
Z+ + Z_ for 1/4 < A' < 1/2, where Z± = tr{e-''^±} 
and Z± = tr{e~'^^±}, respectively. Following the same 
reasoning used in Ref. [l3l it can be shown that Z^ = Z_ 
and = Z_. The impurity entropy of at 
T = has been calculated in Ref. ^ It is In ^/W/2 
for the antisymmetric sector in both H-^- and and 
In \/1K for the symmetric sector in 11 . Together with 
the contribution from . we find that Simy — In (2if ) /2 
for \I2<K < \ and 5„„p = In {2K) for 1/4 < < 1/2. 

The impurity correction to the free energy is given 



ation iujn — )> w + iO"*", where i = y^z and S^'^{t) = 
— {TT{Si{T)Si{0)}) . After performing the unitary trans- 
formation U, we have 



Sy ~ cos 9< sin 



|sin V27rA'<ls(0) 



+ sin 6Sz , 
Sz~ — sin0|sin 



V2^$40) 



Sx + sin 
Sx + sin 



y2^$,(0)]Sy} 
\/2^l>,(0) 5„} 



Near the fixed point, we expect that 5, (r) oc |t|~^' 
at T = as |r| — > -|-oo, and Xi is determined by the 
term in the above equations with the smallest scaling 
dimension,^! Near both the 2CK and the 2CK' fixed 
point, it is the Sx term which determines x^, yielding Xi = 
K for 1/2 < /\ < 1 and 2K for 1/4 < A' < 1/2. Hence, 
we find that at low frequencies Rexi('^,0) « Ci|w|^~^ 
for 1/2 < if < 1 and Ci\uj\'^^-^ for 1/4 < K < 1/2, 
and Rex,(0,T) w CiT^-^ for 1/2 < AT < 1 and 
(j/j^2K-i for 1/4 < if < 1/2 at low temperatures, where 
Cy,Cy (X cos^ 9 and Cz,Cz oc sin^ 9. 

Quantum critical regime. — From the above analysis, 
the difference between the 2CK and the 2CK' phases lies 
at the Hamiltonian for the symmetric sector. Hence, the 
critical theory can be described by the Hamiltonian for 
the symmetric sector with the LL parameter K = 1/2. 
It turns out that this Hamiltonian can be refermionized 



by 6Fi, 



Fo 



jQdTC2{T) +0{X^), where 



Fo is the bulk free energy at the fixed point, C2(t) 

2A 



{T{SH{t)5H{0)})o oc 



with 



being 



sin (ttt / P) 

the expectation value at the fixed point, 6H denotes 
the leading irrelevant operator (LIO), A is the corre- 
sponding coupling constant, and A is the scaling di- 
mension of 5HM Near the 2CK fixed point, the LIO 



is the ^^i^5^sin 



term with scaling di- 



mension K/2 + 1/{2K) for l/VS < if < 1 and 



term with scaling dimen- 



sion 2K for 1/2 < 

in the absence of the Rashba coupling, the LIO is 



as 



K < 1/73. (We notice that Hs^ v dxltpjidxi'i - iplidxA) + m\ii'}ipr{0) 



H.c. 



the ^^^SxCOS Y if"0a(O) term with scaling dimension 

1/if .) On the other hand, near the 2CK' fixed point, the 
LIO is the A2 term with scaling dimension l/(2if ). From 
the above results, we may obtain the temperature depen- 
deuce of the impurity specific heat Cimp = —T-^j^5Fimp 
at low temperatures, yielding Ci^p oc 
l/\/3 < if < 1, r4^-2 for 1/2 < if < 1/V3, Ti/^-2 for 

1/4 < if < 1/3 or 1/3 < A' < 1/2, and T\n[j^^ for 
if = 1/3, where tq the short-time cutoff. 

If we apply a small bias across the two edges, a cur- 
rent will flow from one edge to the other. Since only the 
J2 terms will contribute to this current, the leading tem- 
perature dependence of the conductance G at zero bias 
reflects the renormalization-group (RG) flow of the J2_l 
term. From the above discussions, we see that neither 
the qualitative RG flow of the J2± term or its scaling 
dimensions near the various flxed points are affected by 
the the presence of the Rashba coupling. As a result, the 
temperature dependence of G is identical to that without 
the Rashba coupling. (Sec Fig. 1(b) in Ref. [isl ) 

One way to distinguish the 2CK and 2CK' fixed 
point is to investigate the dynamical structure factor 
of the impurity spin Xii'^^'T), which can be obtained 

from SiiiuJn) = J^dre^'^^'^ sI'^\t) by analytic continu- 



where m = ±dJi±£,a/2 and V//r « e=F»v^(*»±e' 



(8) 
In 

Eq. ([8]), we have set Sx — ±1/2 because [Sx,Hs] = 0. 
Hg is nothing but the Hamiltonian of a one-dimensional 
Fermi liquid backscattered by a J-function like poten- 
tial at X = 0. Therefore, we expect that its thermody- 
namical properties should resemble those of the Fermi 
liquid. In fact, straightforward calculations show that 



is) 
irnp 



(T = 0) and Gl'J^p oc T. Since the combined 



^> - . . 

contribution to Simp from the antisymmetric sector and 
Sx vanishes, we conclude that Simp = at the QCP. 
Moreover, the contribution to Gimp arising from the an- 
tisymmetric sector is given by the LIO, which will give a 
higher power in T . Consequently, the leading tempera- 
ture dependence of Gimp is dominated by the symmetric 
sector, i.e. Gimp oc T. Logarithmic corrections to this 
result arc possible and can be obtained by the one-loop 
RG equation near the critical point if — 1/2. However, 
this is beyond the scope of the present work. 

Concluding Remarks. — To summarize, we have shown 
that although the low-temperature physics of a QD cou- 
pled to two helical edge states is still described by the 
2CK fixed point in the presence of a Rashba coupling for 
1/2 < K < 1, the leading temperature dependence of 
thermodynamical quantities is drastically changed due 
to the new LIO's produced by the Rashba interaction. 
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With stronger Coulomb repulsions between electrons on 
the edges, the system will be driven to a new phases 
described by the 2CK' fixed point, which can be dis- 
tinguished from the 2CK fixed point by examining the 
impurity spin susceptibilities. The two phases are sepa- 
rated by a QCP at K ~ 1/2. Hence, near the QCP, there 
should be a window to observe the quantum critical be- 
havior different from those in the two phases, as we have 
shown. Since the LL parameter K depends sensitively 
on the Rashba interaction strength, the new phase de- 
scribed by the 2CK' fixed point and the quantum phase 
transition into this phase can be electrically controlled. 



Therefore, our results not only reveal a new 2CK fixed 
point that has not been analyzed before, but also serve 
as a useful guide for future experimental investigations 
on this system. 
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